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We study a model of a magnetic resonance force microscope (MRFM) based on the cycHc adiabatic 
inversion technique as a high-resolution tool to detect single electron spins. We investigate the 
quantum dynamics of spin and cantilever in the presence of coupling to an environment. To obtain 
the reduced dynamics of the combined system of spin and cantilever, we use the Fey nman- Vernon 
influence functional and get results valid at any temperature as well as at arbitrary system-bath 
coupling strength. We propose that the MRFM can be used as a quantum measurement device, i.e., 
not only to detect the modulus of the spin but also its direction. 

PACS numbers: 03.65.Yz,07.79.Pk,07.10.Cm 



I. INTRODUCTION 

Magnetic resonance imaging technologies (MRI, NMR, 
ESR) are widely used to characterize physical, chemical, 
and biological samples. What makes them powerful is 
that they are non-destructive and capable to probe the 
three-dimensional structure of the sample^. Recently, 
looking at structures at the molecular or atomic level has 
become important in a number of scientific disciplines. 
Magnetic resonance force microscopes (MRFMs) have 
been developed to bring magnetic resonance imaging 
technologies to such an ultimate resolution. The MRFM 
combines conventional magnetic resonance technology 
with probe microscope technology, e.g., atomic force mi- 
croscopy, to image individual molecules or atoms^. In 
an MRFM, a magnetic particle mounted on a cantilever 
interacts with nuclear or electron spins in the sample via 
the very weak magnetic dipole force. When modulated at 
resonance with the cantilever oscillation frequency, even 
the weak magnetic force induces sufficiently large vibra- 
tions of the cantilever. By probing the resulting vibra- 
tional motion of the cantilever, it is in principle possible 
to detect spins with molecular or atomic resolution. The 
cyclic adiabatic inversion (CAI) technique has been pro- 
posed^ as a promising method to modulate the magnetic 
force. 

The future of the MRFM depends crucially on the de- 
velopment of proper mechanical micro-resonators, e.g., 
cantilevers^. Remarkable progress has been made in this 
direction and the detection of atto-newton or subatto- 
newton scale forces has been achieved already^' ^. Re- 
cently, a nanomechanical flexural resonator at microwave 



frequencies has also been realized^. The development of 
the proper technology to detect nanometer-scale mechan- 
ical motion is also important. Optical interferometry or 
electrical parametric transducers are the most common 
examples^'^'^. In recent work, a single-electron transis- 
tor capacitively coupled to a nanomechanical resonator 
has been used to detect the vibrational motion of the 
resonator even in the quantum regime^. 

The progress in MRFM and related technologies has 
also attracted theoretical interest, especially, the ques- 
tion of single-spin detection using the MRFM. Mozyrsky 
et al.^ studied the relaxation of a spin, treating the can- 
tilever as a classical noise source. Berman et ali^^i^ 
studied a CAI-based MRFM and treated both the spin 
and the cantilever as quantum systems that are subject 
to environmental effects. They addressed two interest- 
ing and important issues: which component is measured 
in an MRFM single-spin measurement and whether the 
two spin states (up and down) lead to distinctively dif- 
ferent cantilever motions. They solved numerically the 
time-dependent Schrodinger equation for the spin-plus- 
cantilever system in the absence of coupling to the en- 
vironment. In the presence of an environment, they 
constructed a generalized master equation in the high- 
temperature limit, and solved it numerically. We note 
that their master equation is based on the Markov ap- 
proximation, and is not in Lindblad form.^^'^^ (the nor- 
malization and the positivity of the density matrix are 
not guaranteed). 

In this paper, we study the measurement of single spins 
with the MRFM based on the CAI technique. The start- 
ing point of our work is closely related to Refs. llOllllL 
In the absence of the coupling to the environment, we 
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FIG. 1: MRFM measurement device. A cantilever carrying 
a magnetic particle is subject to a static magnetic field By 
in the ^-direction, and a time-dependent field B^(t) rotating 
with frequency ojrf in the x-y plane. The cantilever is coupled 
to a sample spin by a magnetic force 77. 



solve the time-dependent Schrodinger equation exactly 
and confirm the numerical results by Berman et ali^^i^ 
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I.e., we 



We use an open quantum system approach 
take the influence of the environment into account by 
coupling a harmonic oscillator bath to the cantilever. To 
calculate the dynamics of the spin during the measure- 
ment process, we take an effective-bath approach, and 
obtain the exact solution for the reduced density matrix 
of the spin. To find the cantilever dynamics, we solve 
the Feynman- Vernon influence functionaUSiii in order to 
obtain the reduced density matrix of the spin plus can- 
tilever system. Both methods are are valid at any tem- 
perature as well as for arbitrary coupling strength (within 
the CAI-scheme). This analytical approach allows us to 
interpret the results in a transparent way and to inves- 
tigate the issue whether the MRFM can be used as a 
quantum measurement device to probe the spin state. 

The paper is organized as follows: In Section ^ we 
first introduce the model and discuss our adiabatic Born- 
Oppenheimer approximation scheme in connection with 
the CAI-technique. In Section Hill we present the exact 
solution of the time-dependent Schrodinger equation for 
the spin-plus-cantilever system without coupling to the 
environment, the results of which will be compared with 
those in the dissipative case in the later sections. In Sec- 
tion llVl we investigate the quantum dissipative dynamics 
of the spin alone using an effective-bath approach. The 
dynamics of the cantilever is investigated in Section IVl 
The physical implications of the solution are analyzed in 
detail and the possibility to use the MRFM as a quantum 
measurement device is discussed in Section IVll Finally, 
in Section rvTll we draw our conclusions. 



II. MODEL 

We consider an MRFM setup based on the cyclic adi- 
abatic inversion (CAI) technique (see Fig. It con- 
sists of a ferromagnetic particle mounted on the tip of 
a cantilever, a strong static magnetic field By in the z- 
direction, and an rf field Bx(t) rotating with frequency 
cjrf in the x-y plane modulated by 0(t): 



COS (cJrft — (/>(t)) 

— sin (cjrft — (pit)) 



(2.1) 



As in usual NMR setups, one puts uo^i = Cz = ^/iB^y, 
where g is the ^-factor of the spin and jj^b the Bohr 
magneton. For later use, we also define e± = qiibB^. 
The "sample" consists of a spin interacting with the fer- 
romagnetic particle via the magnetic force 77 and with the 
static and rf fields. The Hamiltonian of the spin and the 
cantilever is given by 



V<^zZ + ^ + -, (2.2) 



where the a's are Pauli matrices, a± = {ax ± iay)/2, 
and z (pz) is the position (momentum) operator of the 
cantilever. In Eq. (|2.2|) and hereafter we use a unit sys- 
tem such that h = ks = ^0 = ^0 = I, where uj q is the 
natural frequency of the cantilever and Iq = ^JUfrmJ^ 
is the harmonic-oscillator length. It is convenient to 
move to a frame rotating with the rf field by making 
a transformationi& 



(2.3) 



with = exp {-^[^zt — (l){t)]d-z}. The resulting Hamilto- 
nian readsiSiii 



The idea of the CAI-based MRFM is as follows: The 
phase modulation (j){t) of the rf field is assumed to be har- 
monic and causes adiabatic inversions of the spin, which 
in turn exert an oscillating force on the cantilever. At res- 
onance, i.e., if the frequency of the modulation is equal 
to the natural frequency of the cantilever (which is 1 in 
our units). 



4>{t) = ^0 sin(t - (p) , 



(2.5) 



the vibration amplitude of the cantilever can be large 
even for a very small magnetic force rj. 

Equation (|2.4|) describes a spin which couples to a har- 
monic oscillator and is itself subject to a time-dependent 
effective magnetic field gfiB^efi{t) = ^±^x + (p{t)^z^ 
where e^^, Gz are unit vectors in the rotating system. The 
Hamiltonian in Eq. (|2.4|) is not exactly solvable. Here we 
make a plausible approximation based on the following 
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observations. For typical experimental parametersiSiii, 
Beff varies slowly compared with the Rabi oscillation fre- 
quency: |Beff(t)|/|Beff(t)| < e{t) = ^ei + (/)2(t). Ac- 
cording to the adiabatic theorem.^^'^^, the spin part of 
the solution should be determined by the adiabatic evo- 
lution; i.e., the spin "follows adiabatically" the effective 
field Beff(t). It is therefore convenient to choose the ba- 
sis states |x+(^)) and |x-(^)) quantized along the axis 
parallel to Beff(t) (notice that there is no Berry phase 
because the solid angle enclosed by Beff(t) is zero). In 
this basis, the Hamiltonian in Eq. (|2.4|) is recast to 



^{t) = -^^{t)rz - ^tH^^^ + ^^^^^ 

4- 



(2.6) 



where Tx and are the Pauli matrices with respect 
to the frame rotating adiabatically with Beff(t). We 
have suppressed the time arguments to the Pauli matri- 
ces f,{t) = Ix+W) (x+WI - Ix-W) (X- Wl and f,(t) = 
|X+W)(X-WI + |X-W)(X+WI in Eq. This is 

because within the adiabatic approximation, the dy- 
namics of the spin part of the wave function is com- 
pletely governed by the basis states |x±(^)) and the dy- 
namic phases, i.e., |x(^)) = c+e~*-^o ^^'^+(^') + 
c_e~^5o^^'^-^^'^ \X-(t)). We further note that the spin 
dynamics is much faster than the cantilever motion, 
e(^) ^ ex ^ 1- The situation is reminiscent of the Born- 
Oppenheimer approximation^^, where the nuclei interact 
with the average charge density of the electrons which 
move much faster. In our system the nuclei correspond 
to the harmonic oscillator which is interacting with the 
averaged motion of the spin. Therefore one can drop the 
third term in Eq. (|2.6^) . (The deviation of the spin due to 
this term is also negligibly small since ri\ {z{t)) \ <C e(t); 
see below). Using this approximation we finally get the 
following Hamiltonian, which is the basis of the further 
considerations in the paper: 



.rjf{t)f,z^^^^^ (2.7) 



where f{t) = (j){t)/e(t). This form is justified in a more 
rigorous way in Appendix also taking into account the 
influence of the environment (see below). Its validity was 
also confirmed by the numerical simulations in Ref. 0. 

So far we have described a model for an idealized sys- 
tem of spin and cantilever. In reality they are coupled to 
various environments, which lead to decoherence as well 
as damping. In particular, the cantilever is inevitably un- 
der the influence of phonons or other vibrational modes 
which are close in frequency to the single mode in ques- 
tion. The (direct) environmental effects for the spin, e.g., 
hyperfine interaction, spin-lattice relaxation, etc., are rel- 
atively small. Therefore, for simplicity, we assume a sim- 
ple Ohmic bath of oscillator s^^i^^i^^i^^ directly coupled to 



the cantilever but not to the spin. Then the total Hamil- 
tonian for the spin and the cantilever plus the oscillator 
bath is given by 



^otal(^) = ^(t) 

CO 



2m/e 



rrikujl 



Ck 



(2.8) 



All the relevant features of the Ohmic bath are charac- 
terized by the spectral density 



= aujQ{l — uj/ujc) • 



(2.9) 



where a is a dimensionless parameter characterizing the 
coupling between the system and the environment and 
coc is the cut-off frequency. The spin dynamics and the 
probability distribution of the cantilever will not depend 
on the cut-off. 

We describe the system of spin plus cantilever in terms 
of the reduced density matrix p{t) = tr^ptot(^) by tracing 
out the bath. In the realistic typical experimental situ- 
ation, the cantilever always remains in contact with the 
environment. Thus, the cantilever and bath are not in 
a product state at the beginning of the experiment. For 
the calculation with the influence functional, we can take 
this fact into account, assuming that the cantilever and 
the bath were in a factorized state at a time t = to. In 
the limit to ^ — oo we get then the realistic initial state 
for the cantilever at the time t = 0. If we would start 
with a factorized state between cantilever and bath, the 
solution would be very sensitive to the initial condition 
of the cantilever; see Section HVl 

Furthermore, it is assumed that the interaction be- 
tween the spin and the cantilever is turned on at t = 0, 
i.e., f{t) = for t < 0. The measurement happens at 
times t > 0. The initial state p(0) of the density matrix 
is a product state, 

p(0)=/5(^)(0)p(^)(0), (2.10) 

where p'^^^ is the density matrix for the spin only and 
p^^^ describes the cantilever in thermal equilibrium with 
the bath. From the CAI scheme and from the associated 
adiabatic approximation discussed above it then follows 
that the density matrix at times t > has the form 

{s,z\mW,z')=pi^}{0)p^^;{z,z',t). (2.11) 

Namely, the dynamics of the density matrix p{t) is com- 
pletely determined by the spin-dependent cantilever part 
p[^){z,z',t). 

Here the spin-dependent cantilever part should not be 
confused with the density matrix for the cantilever only, 
which is given by 



p(^\z,z',t) = Y,{s,z\p{t) 



5, Z 



S = ± 



pV+(0)pS(^, z\ t) + pL^2(0)pL^i(z, z', t) . (2.12) 
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Analogously, the density matrix for the spin only at time 
t > is given by 



pi'J{t)=plViO) I dzpl'^/iz^z^t). (2.13) 



There are several ways to prepare the spin in a partic- 
ular state^^, and we will assume a general state p[^}{0). 



III. THE COHERENT SOLUTION WITHOUT 
BATH 

Before we investigate the full Hamiltonian in Eq. (|2.8|) , 
it will be instructive to first consider the problem with- 
out bath, Eq. (|2.7|) . The time-dependent Hamiltonian in 
Eq. (|2.7|) can be solved exactly for arbitrary functions e{t) 
and f{t) of t (of course, the variation of e(t) and f{t) in 
time should be sufficiently slow so that the Hamiltonian 
Eq. (|2.7|) is meaningful). 

One can show that the time-evolution operator 

^(t2,ti) = fexp -i J^^ dt'je{t') (f is the time- 
ordering operator) is given by 



^(t2,ti) =exp 



where 



X ^{r,at2))%{t2 - ti)#(f,^(ti)) , (3.1) 



^(t)^zrj^l^dt'e-^^'-''^f{t'), (3.2) 



^o(^) = exp {-itd^d) 



(3.3) 



a = {z + ipz)/ V^, and ^{(,) is a displacement operator^'* 
defined for a complex number ^ by 



^(0 = exp(Cat _ . 



(3.4) 



The coefficient 0(^1,^2) in Eq. is a real function of 

ti and t2 (one does not need an explicit expression of it 
because it drops out of the following calculations). 

To illustrate the dynamics created by the time- 
evolution operator in Eq. (|OJ, let us discuss an example. 
Suppose that we start at time t = with the cantilever 
in a coherent state 



: exp 



(3.5) 



and with the spin in a linear superposition (with ampli- 
tudes c+ and c_) 

|x(0)) = c+|x+(0))+c_|x-(0)). (3.6) 

The total wave function at t = is given by 

|f(^,0))=V(^,0)|x(0)) , (3.7) 



and at later time t > 0, by 

\<i>{z,t)) = c+i;+{z,t)\x+{t))+c.i^-{z,t)\x-{t)) . (3.8) 

The cantilever wave function in Eq. (|3.8|) for each spin 
component is given by 



i^±{z^t) = ^exp 



ic{t 



,0)±i [ dt' e{t') 
Jo 



xexp<,--z'^V2e^{t)z-[Ree^{t)]'\ , (3.9) 



where 



(3.10) 



Therefore, the average position of the cantilever is 
{z{t))^ = >/2Re^^(t) for spin s = ±, respectively, 
whereas the average momentum is given by {pz{t))_^ = 

V^Im^^(t). Here it is interesting to note (in compar- 
ison with the results below) that exactly at resonance 
[see Eq. (IT^ ]. |^(t)| in Eq. [and hence \^±{t)\ in 

Eq. (|3.10|) ] contains a term which linearly increases with 
time t. In other words, the oscillation amplitude of 
the cantilever gets indefinitely larger and larger as time 
passes. This is not surprising since we are driving an ideal 
oscillator at the resonance frequency, and in fact this is 
what allows the MRFM to detect ultra-small forces. In 
reality, the cantilever is subject to various environmen- 
tal effects and the oscillation amplitude is bounded from 
above (i.e., the Q-factor is finite). This is the case that 
we will study below. 



IV. DYNAMICS OF THE SPIN 

Now we take into account the influence of the bath. In 
this section, we ffist analyze the dynamics of the spin. 
The dynamics of the cantilever will be discussed in the 
following section. When we are interested in the dynam- 
ics of the spin alone, we can regard the cantilever as a 
part of the environment. In fact, Garg et al?^ (see also 
Refs. 12602711281 ) showed that the problem is equivalent to 
a spin coupled linearly to an oscillator bath: 



k 

-^^cokblbk. (4.1) 

k 

The distribution of the oscillator frequencies cok and the 
coupling constant gk are now characterized by a non- 
Ohmic spectral density 



^ — ^ TT [00^ — 1 



(C^2 _ 1)2 + (^^)2 



(4.2) 



To investigate the spin dynamics, we write the reduced 
density matrix of the spin 

p^'Ht) = tVB^tot{t)ptot{0)^tlit) (4.3) 
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in terms of the time-evolution operator "^^tot {t) associated 
with J^tot{t) in Eq. (^^. In analogy to Eq. (|OJ, the 
time-evolution operator is given by 



^totW =exp 



k 

(4.4) 



where 



[ dt't 
Jo 



i{t-t')ujk 



fit'), 



(4.5) 



and ^ is now the displacement operator for the k-th 
mode of the bath, i.e., d should be replaced by bk in 
Eq. 

For the initial state ptot(O), we assume [see Eqs. (|2.1()|) ] 



Ptot(o)=p(^)(o)n 



-f3ujkblbk 



(4.6) 



Then the density matrix for the spin is given by 



As-s') 



f dt'eit') 
Jo 



xn(^^(^'^'oW)^(^?'oW))fe' (4-7) 

k 

where (•••)/. is the average with respect to the k-th os- 
cillator in the bath. 

Equation (|4.7|) shows that the diagonal elements of the 
density matrix {s = s^) are constant in time 



(0). 



(4.8) 



In other words, there is no spin relaxation and the spin 
dynamics is pure dephasing because there are no trans- 
verse fields. This is consistent with the adiabatic approx- 
imation we made at the beginning. 

On the other hand, the off-diagonal elements {s ^ s') 
are expected to vanish rapidly with time. This can be 
seen from [see Eq. (|4.7|) ] 



T{t)^i I dt'e{t') 





pV^2(t)=pi^2(0)exp 
where 

r(t)^2^|a(i)pcoth(^), 

k 

or in terms of the spectral density function 



(4.9) 



(4.10) 



dt'e'"^'' fit') 



T{t) = 2r/2 j Jeff (w) coth 

(4.11) 

Figure 121 shows {t)\ evaluated using Eqs. (|4.9|) and 

(|4.11|) . To compare our results with those of Berman et 
a/iMl who assumed an initial product state of cantilever 




FIG. 2: Main plot: (t)| for different temperatures T — 

0, 1, 2, 5, 10, 100, for 0o = 1000, e± = 400, ry = 0.3, a = 0.006, 
and ix}c — 1000. The initial condition for cantilever and bath 
is the thermal equilibrium state. Inset: same quantity for 
an initial product state of cantilever and bath. Initially, the 
cantilever wave function is a Gaussian with width a — \f2. 
In both cases, pif?(0) = 1/2 for s, = ±. 



and bath, the inset of Fig. [2 shows for a Gaussian 

initial state of the cantilever. (To obtain these results 
we evaluate the path-integral formulas in Appendix El 
with to = instead of taking the limit to ^ — oo ). If 
we compare the main part of Fig. [21 with the inset the 
strong dependence on the initial conditions is evident. 
The slower decay and the more pronounced oscillations 
shown in the inset are a consequence of the oscillatory re- 
laxation of the cantilever to its thermal equilibrium state 
if one starts with an initial product state of cantilever 
and bath. On increasing the coupling a, the oscillatory 
behavior becomes less visible since the cantilever relaxes 
immediately to its thermal state. 



V. DYNAMICS OF THE CANTILEVER 

In Section iTTTl we described the driven dynamics of the 
otherwise isolated system of spin and cantilever deter- 
mined by the Hamiltonian Eq. (|2.7|) . In this section, we 
now take into account the inffuence of the environment 
starting from the Hamiltonian Eq. (|2.8|) . The reduced 
dynamics is obtained analytically with the Feynman- 
Vernon inffuence functional^Siii for arbitrary coupling 
strength a to the bath and for arbitrary temperature 
T. The advantage of this method as compared Ref. O 
is that no master equation is used and that there is no 
restriction on the number of basis functions used to nu- 
merically integrate the problem. 

The reduced dynamics of the cantilever obtained with 
the inffuence functional is given by: 

pfj{zs,z),t)= (5.1) 
= y dzidzlJs.,'{zf, z'f, t; Zi, z\, to)p'^J{zi, z\, to) , 
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where the influence functional is 



s^s' = ±, and the action is deflned by: 
S,4z,z']=S''Az]-S''Az'] 



(5.2) 
(5.3) 



[ dT[z{T)-z'{T)][z{T)+i'{T)] 
Jta 



I to 

rt ft 



+ U dr [ dT'[ziT)-z'iT)]KiT-T')[ziT')-z'iT')] 
^ Jto J to 

This form of the action is only valid for an Ohmic bathSi. 
Furthermore, K(t) is the real part of the bath correlation 
function 



K{t)= Re(l(T)X(0)), 
where X{t) = "^^CkXkit). Finally, 



(5.4) 



= dT[-z\T) - -z\t) + r,sf{T)z{T) + -e{T)s] 

(5.5) 

is the bare action without oscillator bath. 

The action can be simplifled further by introducing 
relative coordinates deflned by R = {z -\- z')/2 and r = 
z — z'. The action is then found to be 

Sss' [R. r] = Si. [R, r]-a [ drR{r)r{r) 



dr 



to J to 



[ dr'r{r)K{ 

Jta 



r-ry{r'), (5.6) 



with 



= I dT\R{T)f{T)-R{T)r{T)+r]f{T)R{T){s-s') 

Jto ^ 

+ \vf{r)r{T){s + s') + ie(r)(s - s')} . (5.7) 

In the next step, the action is expanded around the 
classical path. The classical equations of motion can be 
found by minimizing this action and read 

R{r) + aR{r) + R{r) = F^(r) , (5.8) 



f(r) — ar{r) + r(r) = Fr{r) , 



(5.9) 



FR{r) = lr]f{T){s^s')^i f dr' K{t -T')r{T') , (5.10) 

^ Jta 



with classical solutions Rci{t)^ '^cii^)^ respectively. Note 
that the solutions are complex^, and the dependence 
on 5, 5^ of all these quantities has been suppressed. The 
classical solutions, which are given in Appendix El ^ire 
linear in the boundary values i?/,r/,i?i and r^. There- 
fore, Sss' [Rchrd] is a bilinear form in these variables. We 
obtain 

Jss'{Rf,rf,t;Ri,ri,to) = ex.p (^iSss'[Rci,rci]^ , 

(5.12) 

where all the contributions from the fluctuations 
around the classical path are contained in the time- 
dependent, but spin-independent normalization constant 
Af{t) , which can be obtained from the normalization con- 
dition 



/CO 
dRfp,,{Rf,rf = 0,t) = l. 



(5.13) 



The Gaussian form of the expressions leads to a flnal 
reduced density matrix of Gaussian form if the initial 
density matrix is Gaussian, which is true for a coher- 
ent state. Therefore we deal with Gaussian wave packets 
also in the dissipative case. The explicit formulas are dis- 
cussed in detail in AppendixEl where the solution for the 
reduced dynamics is obtained starting from a Gaussian 
wave packet at time to. We then take the limit to ^ —oc 
such that the information about the initial state is lost 
at time t = 0. 

We will now give analytical expressions of the density 
matrix for the diagonal and off-diagonal elements with 
respect to the spin degree of freedom. Let us flrst discuss 
the result for s = s^: 



pfs\R.r,t) 



^^v{-i^[R-Xs{t)f 

R 

^irxs{t)}, (5.14) 



2tT(Jr ^ ^ 2cr2 



where the flnal coordinates have been replaced hy R = 
Rf and r = Vf. The widths of the Gaussian peaks are 
independent of the spin. The width in the i?-direction is 
given by 

cr| = ^ duJ^^iuj) coth (^^^ . (5.15) 

cFji increases with temperature. This is because the can- 
tilever position suffers more thermal fluctuations. The 
width in r-direction is found to be 



~2 ^ J diJu;'^Jefi{u;)coth(^^^. (5. 



16) 



Fr{r)=r]f{r){s-s'). 



(5.11) 



Note that as is well-known the momentum width diverges 
with the cut-off frequency coc which was deflned after 
Eq. (|2.9|) . That is why we kept the dependence on the 
cut-off in this integral. The spin dynamics and the prob- 
ability distribution of the cantilever will not depend on 
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the cut-off. In contrast to cfr^ cr^ decreases with temper- 
ature; this is natural since the cantilever gets closer to a 
classical oscillator as temperature goes up. The temper- 
ature behavior of these two integrals can be read off in 
the limit of small a <C 1, viz., 



4 



1 



-coth . 
2 \2TJ 



(5.17) 



The Gaussian wave packets are moving according to 

rt sin [ujR{t — t') ) 
Xs{t)=r]s dt'e-^^'-''^ ^ -f(t'), (5.18) 

which d epends on t he spin s. The oscillator frequency 
ujR = Y^l — (q^/2)^ is renormalized due to the coupling 
to the bath. Furthermore, Xs{t) is the solution of the co- 
ordinate of a classical dissipative driven harmonic oscil- 
lator with a spin-dependent driving force r]sf{t) starting 
from the initial conditions Xs{0) = and Xs(0) = 0. So 
the result becomes very clear, because the classical so- 
lution is well-known to be an oscillating function, which 
goes trough a transient regime and for t 1/a the am- 
plitude of the oscillation saturates at a finite value. The 
oscillation is periodic (but not necessarily sinusoidal) in 
time with unit period (Tq = 27r/uJo). Consequently, for 
t 1/a the density matrix will show a generic steady- 
state behavior independent of the details of the initial 
preparation of the system. 

The density matrix pi^\R^r^t) behaves quite differ- 
ently with respect to the coordinates R and r. As a 
function of i?, piJ^ (i?, r, t) is a Gaussian distribution with 
standard deviation an and average {R{t)) = Xs{t). On 
the other hand, Xs{t) is the velocity of a classical oscilla- 
tor [see above] , it shows oscillatory behavior in t and r su- 
perimposed on the Gaussian envelope with width a^; see 
Figs. IHHni Thus, the off-diagonal elements pi^\z^z'^t) 
{z z') exhibit an oscillating behavior in t. However, 
this should not be confused with a coherent oscillation, 
which is not expected in this long-time limit. The oscil- 
lation is a consequence of the external driving f(t) [i.e., 
frequency modulation ^(t)]. The diagonal elements (both 

in s and z) pis (^, ^, t) do not show such an oscillation. 

The behavior of Xs{t) can be illustrated by approxi- 
mating f{t) by its primary oscillation amplitude: 



f{t) ^ /o sin(t) + (higher harmonics) , (5.19) 



where 



/o 



^ (g) [Ei-<Pl/el) - Ki-c^l/el)] . (5.20) 



Here, K{x) and E{x) are the complete elliptic integrals 
of the first and second kind^^. One obtains 



cos(t) 



]) 



-(higher harmonics) . 



2ujR 
(5.21) 



This solution shows the main features of the spin- 
dependent separation Xs{t)^ namely the transient be- 
havior and the steady-state oscillations: Xs{t) ^ 
—r]sfocos{t)/a. It is interesting to notice that the av- 
erage cantilever motions are exactly in opposite phases 
(shift by tt) for spin up {s = +1) and down {s = 
— 1). This was also concluded from the numerical 
simulation presented Refs. 10 11. Thus, the MRFM 
can be used as a quantum measurement device, i.e., 
to detect the state of the spin; see below. There- 
fore, if we start initially with the two spin components 
populated, p^^^(O), pL^l(O) > 0, then p^^\R,r,t) = 

p^_^l({))p^^l{R,r,t) p^fl{0)p^_^l{R,r,t) will show two 
peaks moving in opposite directions as time goes on; see 
discussions above and Figs. (SHHl It should be stressed 
that to separate the two peaks with sufficient resolution, 
the widths of the peaks, Eq. (|5.15^) . should not be larger 
than the maximum separation, rif^/a] see Eq. (|5.2H) . 
Clearly, this criterion restricts the maximum operation 
temperature of the device. Figures IHHSl show the typical 
behavior of the density matrix p^^^ (i?, r, t) of the can- 

tilever for p\J(Q) = 1/2 for = ± as initial state. As 
the coupling to the environment a increases, the distance 
between the peaks shrinks and they are harder to distin- 
guish; see Figs, inland El A similar behavior is observed 
as the temperature increases with a fixed; see Figs.jSland 



Now we turn to the off-diagonal elements s = —s^: 
1 



27TaR 



2a? 



where 



^rCs{t)-T{t)^i dt'e{t')}, (5.22) 



,{t) = 2r]sJ^ dcjJeff(cj)coth (5.23) 



X / dt' f{t')co^{u{t-t')) 
Jo 



and 



Cs{t) = 2r]s du;u;Jei^{u;)coth (5.24) 



/ dt'f{t') sin {ij{t-t')) . 
Jo 



(C) 

In r-direction, pl _g{R,r,t) has a Gaussian shape cen- 
tered at (s{t)/(j'^ with width a^. In i?-direction, it is an 
oscillatory function imposed on a Gaussian envelope with 
width aR. Overall, the function pi -^(i^, r, t) decays with 

t in the same way as shown in Fig. 13 i.e., p[^) {R^ r, t) for 
s ^ s' can be observed only in the transient regime. The 
decay is described by the function r(t); see Eq. (|4.1H) . 
Note that a trace over the cantilever dynamics leads us 
back to the results obtained in a much simpler way in 
Section ITVl 
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FIG. 3: \p^^\R,r,t)\ for a time series in the steady-state 
regime starting at time tm at which the two peaks are not 
separated, e.g., tm = 988 The units have been chosen such 
that both the natural frequency uo of the cantilever and its 
harmonic oscillator length are equal to 1. To = 27r/a;o, « = 
0.006, T = 100; the other parameters are as in the caption of 
Fig. 121 The interference fringes are due to the driving. 



VI. MRFM AS A QUANTUM MEASUREMENT 
DEVICE 

One of the conclusions of the analysis presented here 
is that the cantilever oscillates with the same amplitude 
for both initial spin states (up and down). Probing the 
amplitude of the cantilever vibration can only tell the 
absolute value of the spin in the direction of Beff(O), but 
not its sign. However, the oscillations for the initial spin 
up and down states are completely out of phase (phase 
difference of tt) ; see Section This fact was also no- 
ticed by Berman et al}^ in their numerical simulations. 
Hence, there is the possibility to use the MRFM as a 
quantum measurement device, i.e., to detect the direc- 
tion of the spin with the MRFM by probing the (dis- 
crete) relative phases of the cantilever oscillations. In 
the quantum theory of measurement, this falls into the 
category of the indirect quantum measurement scheme^^. 
In such a scheme, the quantum object supposed to be 
measured is coupled to another quantum system, the so- 
called quantum probe. The classical measurement de- 
vice then detects the quantum probe instead of probing 
directly the quantum object. In our case, the quantum 



FIG. 4: \p^^\R,r,t)\ for a time series in the steady-state 
regime for a = 0.012 and T = 100. 



object is the spin, the quantum probe corresponds to 
the cantilever, and the classical measurement device can 
be, e.g., the fiber-optical interferometer. A conceivable 
scheme to measure the relative phases of the cantilever 
oscillations is to use a reference spin which is prepared 
in a definite known state, for example, by applying a 
strong magnetic field in a desired direction. The two sig- 
nals from the reference spin and the spin in an unknown 
state are superposed to determine the relative phase of 
the unknown spin. 



VII. CONCLUSIONS 

We have studied the CAI-based MRFM as a high- 
resolution tool to detect single spins. The quantum dy- 
namics of the spin-plus-cantilever system was analyzed 
in terms of the reduced density matrices, p'^^^t) (for the 
spin) and p^^\t) (for the cantilever), in the presence of 
coupling to the environment. Using an effective bath 
model, we were able to determine the dynamics of the 
spin during the measurement process. Our results re- 
main valid at all temperatures as long as the adiabatic 
approximation is satisfied. We have evaluated the in- 
fluence functional for the combined system of spin and 
cantilever to obtain the quantum dissipative dynamics of 
the cantilever. These results are valid for all tempera- 
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FIG. 5: \p^^\R,r,t)\ for a time series in the steady-state 
regime for a = 0.012 and T = 50. 



FIG. 6: \p 

regime for a 



^^^-^ (i?, r, t) I for a time series in the steady-state 
0.012 and T = 10. 



tures and coupling strengths. Finally, we have proposed 
that the MRFM can be used as a quantum measurement 
device, i.e., not only to detect the absolute value of the 
spin but also to detect its direction. 

The dissipative dynamics of an open quantum system 
is sensitive to the low-frequency behavior of the spec- 
tral density of the environment. While the Ohmic model 
Eq. (|2.9|) is a plausible model, it will be worthwhile to 
identify the sources of the environmental fluctuations and 
construct a physical model of the environment starting 
from a more microscopic theory of the cantilever. 
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APPENDIX A: ESTIMATION OF THE 
SPIN-FLIP RATE 

The cyclic adiabatic inversion scheme implies two ba- 
sic assumptions: (i) The variation of the external driving 
(J){t) is slow enough to allow for an adiabatic approxima- 
tion^^, i.e., <C e^. (ii) The time scales of the spin 

dynamics and the cantilever dynamics are well separated 
{e±_ 1) such that the Born-Oppenheimer approxima- 
tion is justified. Yet, the finite rates of change in the ex- 
ternal driving and the cantilever position will induce spin 
flips. The spin-flip rate can be estimated by the Landau- 
Zener transition (adiabatic transition) rate-— 
For this purpose, we rewrite Eqs. (|2.4|) and (|2.8|) in the 
form 



1 ^ 



(Al) 



where F{t) = (j){t) + 2r]{z{t)). The back-action of the 
cantilever has been accounted for by its time-dependent 
average position, and the contribution from it will be esti- 
mated below in a self-consistent way based on the results 
in Section 13 The probability that the spin flips against 
the effective magnetic field Beff(t) during one period (i.e., 
27r/uJo) is then given by 



Plz ^exp( ^) 



(A2) 



10 



where we have taken v = max 
worst case. 

It fohows from Eqs. (|2.5|) and (I5.21|) that 



V < max 



- 2r] max 



dt 



m) 



to estimate the Note that the back-action of the cantilever is stronger 
for larger Q-factors of the cantilever (Q — l/<^) since the 
maximum velocity of the cantilever increases with the 
Q-factor. 

h + 2^fo. (A3) 

a 



Therefore, assuming typical values for the parameters, 
^0 ~ 1000, e±_ ^ 400, T] ^ I, and a ^ 0.001, we have 
fo^l and 



Plz < exp -TT 



(/)o + 2r]^fo/a 



APPENDIX B: PATH-INTEGRAL FORMULAS 
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-70 



(A4) 



In this appendix we will fill in some of the details left out in Section It is convenient to define 7 = a/2 as the 
friction constant. The classical solutions to Eqs. (|5.8^) . (|5.9^) are given by 



rci{r) 



sin {uR{t -to)) 



[r, sin {unit - r))e^(^-^°) + [vf - r^(t)] sin {u;r{t - to))e^^--'^] + r^(r) , (Bl) 



Rci{r) = —^—^\r^ sin {conit - r))e-^^^-''^ + [Rf - Rp{t)] sin {u;r{t - to))e-^^^-'A + Rp{r) , (B2) 
sm [UR[t - to)) I ' ' ^ ' ) 



where 



rp{T) = r dT'Gr{T-T')Fr{T')., 
•Jto 



(B3) 



R^{t)= r dT'GR{T-T')FR(T'). 



and the Green's functions are defined by 



sin(cji^r) 



(B4) 



(B5) 



Gr{r) =e(r)e' 



sin(cc;i^r) 



The infiuence functional for = 5 is found to be 



(B6) 



JssiRf,rf,t;Ri,ri,to) 

\jm 

2n 



(B7) 



exp [i[Kf{t)Rfrf + Ki{t)Riri - L{t)Rirf - N{t)Rfn + ai{t)n + a/(t)r/] - A{t)r} - B{t)rfn - C{t)rf) , 



where the functions appearing in the influence functional are all real and defined by 

Kf{t) = u!R cot {ujR{t - to)) - 7 , 



(B8) 



Ki{t) = LOR cot {ujR{t - to)) + 7 : 



(B9) 



m 



UlRB 



-id-to) 



sin {ujR{t -to)) 



(BIO) 
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N{t) = . ^ (Bll) 
sm [ujR{t -to)) 



^W = o-^^^^^ / / sin(u;^(r-to))i^(r-rOsin(^^(r'-to))e^^^+^\ (B12) 

2sm [ujR{t-to)) Jto Jto 

-7(t+to) ft ft 

= . 2 / u . ^^ / / (^^(^ - ^))^(^ - (^^(^' - ^o))e^^^+" ^ , (B13) 

sm [u;R{t - to)) Jto Jto 



I g-27to rt rt 

2sin^ {ujR{t-to)) Jto Jto 



a/(t) = x(t)-i^/(tMt), (B15) 

ai(t) = N(t)x(t), (B16) 

x(r) = rjs r dr'GRir - T')f{T') , (B17) 

Jto 



Jtn 



(B18) 



In all of these expressions, the dependence on to has been suppressed. 

Let us now discuss the solution for the density matrix. At time t = to we start in a product state between cantilever 
and bath. The cantilever density matrix is assumed to be a Gaussian wave packet with width cr at t = to, 



27rcr 

One could start from a more general initial state, but we will later take the limit to — oo, such that all the 
information on the initial state is lost completely at time t = 0. The experiment starts at time t = by switching on 
the magnetic field. At this time the cantilever has interacted with the bath for a very long time and is in equilibrium 
with the bath, i.e., not any more in a product state. 

(c) 

The general solution for the diagonal elements of pl^/ starting from this initial condition is 



X exp ■ 



27r /D(i) 

2 



r}\^- A{t) + [2B^{t) - 8A{t)C{t) - -^^^ - 4(^A{t)Kf{t) + L{t)[B{t)Ki{t) + C{t)L{t)]y^ 
+irf |^a/(t) -4[ai(t)B(i) -2a/(t)C(i)]a2 +4ifi(i)[a/(i)ifi(t) +ai(t)L(i)](T4j 
+iRfrf [Kfit) + ^[2C{t)Kf{t) + Bimit)^ + 4Ki{t)[Kf{t)Ki{t) - L{t)N{t)]aA 



-2[ai{t)-N{t)Rf]''a^ 



lD{t) , (B20) 



12 



where 

D{t) = 1 + SC{t)(j^ + ^Kf{t)cT^ . (B21) 

In the hmit to ^ — oo we obtain the final result presented in Eq. (|5.14|) . 
The infiuence functional for s' = —s is found to be given by 

Js,-s{Rf,rf,t;Ri,n,to) = (B22) 



where 



2tt 
X exp 



exp (i[Kj{t)Rfrj + Ki{t)Riri - L{t)Rirj - N{t)Rjri + Aj{t)Rj + Ai{t)Ri + ^ dTe(T)]) 
_ A{t)r) - B{t)rfn - C{t)rl + h{t)n + bf{t)rf + b{t)) , 



Af(t)=y{t)-Ki{t)y{t), (B23) 



Ai{t) = L{t)y{t) , (B24) 



bfit) = 2A{t)y{t) - f dr f dr'y{r')K{r - r') ""^"^^: uK^^^^^ ' (^^5) 
J to J to sm[ijR{t -to)) 



J to J to sm - tojj 



b(t) = -A{t)y\t)+y{t) fdr f dr'y{r')K{r - r' f' ^"^^^ ]f'~T^' 
J to J to sm[ujR(t-to)) 



I f dr f dT'y{T)K{T - T')y{T') , (B27) 

^ Jtc^ J to 



'to to 



y{T) = 27?s r dr'GriT - T')f{T') , (B28) 

Jto 



y{T) = 2tJS f dr'drGriT - T')f{T') . (B29) 

Jto 
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This leads to the fohowing general expression for the off-diagonal elements of p^^J'- 
pl/,,{Rf,rf,t) 



X exp ■ 



2n /D(i) 

A{t) + [2B\t) - 8A{t)Cit) - - 4:(^A{t)Kf{t) + L{t)[B{t)Ki{t) + C(t)L(t)]y^ 



\bf{t) + [-^B{t)bi{t) + 8bf{t)C{t) + Ai{t)L{t)y 
+4[Aiit)Bit)Kiit) + bfit)Kfit) + 2Ai{t)C{t)L{t) + hi{t)Ki{t)L{t)]a^ 
+iRfrf (^Kfit) + 4[2Cit)Kfit) + B{t)N{t)]a^ + AKi{t)[Kf{t)Ki{t) - L{t)N{t)]a^ 
+iRf\ Af{t) +SAf{t)C{t)a'^ + AKi{t)[Af{t)Ki{t) + Ai{t)N{t)]a^\ 



+2[6i(t) - iN{t)Rffa'^ 

Am 



a' - AAi{t)[Ai{t)C{t) + bi{t)Ki{t)]a'' 



/D{t) 



-I dTe{T) + b{t) j. . (B30) 



The reduced dynamics of the spin alone is found by tracing out the cantilever coordinates. The result is 

AS) _ (S) ^^2.^^ C{t) AfjtMt) Ajjt) 

- P.,-.(0)exp( Af{t)^^^^^^ ^^^^ g^2^2(^) (^31) 

[Af{t)K,{t) + A,{t)N{t)]^ + i f dre{T)) 



^pif!,(0)exp(-r(t)+z f dreir)). 



The decay rate r(t); see Eq. (|4.11|) . can be obtained in the limit to ^ — oo after a straightforward but tedious 
calculation. In the same limit, we get the result for the den sity matrix presented in Eq. (|5.22|). 
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